No evidence for anomalously low variance circles on the sky 



O 

(N 

5h . 

< 



u: 
^' : 

Or 
i < 

O 
o3 



> 
in 
o 
m 

(N 

o 



X 



Adam Moss,* Douglas Scott,''' and James P. Zibin* 
Department of Physics & Astronomy 
University of British Columbia, 
Vancouver, BC, V6T 1Z1 Canada 
(Dated: April 8, 2011) 

In a recent paper, Gurzadyan & Penrose claim to have found directions on the sky centred on 
which are circles of anomalously low variance in the cosmic microwave background (CMB). These 
features are presented as evidence for a particular picture of the very early Universe. We attempted 
to repeat the analysis of these authors, and we can indeed confirm that such variations do exist in 
the temperature variance for annuli around points in the data. However, we find that this variation 
is entirely expected in a sky which contains the usual CMB anisotropics. In other words, properly 
simulated Gaussian CMB data contain just the sorts of variations claimed. Gurzadyan & Penrose 
have not found evidence for pre-Big Bang phenomena, but have simply re-discovered that the CMB 
contains structure. 



A recent preprint by Gurzadyan & Penrose [1] (here- 
after GP10) makes what must be considered by any stan- 
dard to be an extraordinary claim - that the cosmic 
microwave background (CMB) contains statistically un- 
usual circles which are signatures of events which hap- 
pened before the hot, dense stage in the usual Big Bang 
theory. In particular, GP10 claim to have found direc- 
tions on the sky centred on which are circles of anoma- 
lously low variance in the CMB. If true, this would have 
profound implications for our understanding of the early 
Universe, since such features are not expected in the stan- 
dard picture of primordial structure, based on Gaussian 
random fields. Such a result would also be very surpris- 
ing, given the intense scrutiny the CMB has seen in the 
search for signs of non-Gaussianity. However, we will 
show that the evidence claimed in GP10 fails to be ex- 
traordinary; indeed, such low- variance circles are expected 
to arise, as properly simulated data show. 

GP10 studied maps from the Wilkinson Microwave 
Anisotropy Probe ( WMAP) 7- year data release [2]. They 
specifically looked at circles around particular positions, 
calculating the variance within annuli of 0.5° width, and 
plotting this variance as a function of radius. They found 
some positions where there appeared to be concentric 
rings of suppressed variance. In addition, figure 2 in 
GP10 implies that the scatter in the annulus variances is 
much larger in some particular locations than it is for a 
simulated CMB sky. 

We first tried to reproduce the plots of the variance in 
annuli presented for the two selected positions in GP10. 
We found that we could reproduce very similar results 
(see Fig. 1), and that these are present in both W-band 
and V-band maps, as well as in the so-called ILC map (see 
[3]), which is largely free of foreground signals. Hence it 
is clear that the particular map used is not the issue, and 
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FIG. 1: Standard deviation in temperature measured within 
circular annuli bounded by radius R and R + dR, with dR = 
0.5°, as a function of R for the WMAP 7- year W-band data. 
The red and blue curves are for the two specific directions 
selected by GP10 (red corresponds to galactic coordinates I = 
105°, b = 37°, and blue to I = 252°, b = -31°), and agree 
well with the results presented there. The dashed curves are 
for 500 randomly chosen centres (only 100 are shown). The 
green curve was chosen from the random centres on the basis 
that it has a prominent low-variance annulus at R = 2.5° as 
well as another near R = 6° . Other curves having annuli with 
high variance could just as easily be highlighted. Note also 
that the vertical axis on this plot does not start at zero, so the 
variations in variance are only at the ~ 10% level, and hence 
are quite subtle in the maps themselves. 



it is also unlikely that foregrounds are a major part of the 
effect. Moreover, we found that strong variation among 
the annuli could be seen all over the map. In other words, 
the particular positions selected in GP10 did not seem 
very unusual. Additionally we could find rings of high 
variance just as easily as those of low variance (examples 
can be seen in Fig. 1); this is just what would be expected 
if both low and high variance rings are simply random 
excursions. 
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For all of the results presented here we used the masked 
CMB maps. In each case we choose 500 random direc- 
tions at Galactic latitudes of |6| > 30°. Also, because 
of the WMAP scanning process, the number of samples 
going into individual pixels can vary, and so the uncer- 
tainty in the temperature varies across the sky. Therefore 
our results include proper noise weighting and removal of 
bias. As a check, we also performed the calculations us- 
ing the unmasked sky, and not taking into account noise 
weighting. In this case, our results for the two directions 
chosen in GP10 matched those of GP10 almost perfectly. 
However, ignoring these details did not have a large effect 
on our results (despite one of the directions suggested in 
GP10 being very close to the north ecliptic pole). 

We next turned our attention to simulations. We 
performed the same variance calculations as described 
above, but applied to a combination of simulated maps 
for all W-band feedhorns over the full 7 years (as pro- 
vided on the LAMBDA website [3]); this should represent 
a realistic simulation of the CMB signal plus the WMAP 
noise. Our results are presented in Fig. 2. Comparing 
with Fig. 1, it is clear that the expectation of the annu- 
lus variances, together with their degree of fluctuation, 
are indistinguishable between the simulated data and the 
real data. Since the WMAP data are dominated by the 
CMB anisotropics rather than by instrumental noise, it 
is not surprising that we found qualitatively the same re- 
sults using simpler simulations which were based purely 
on the estimated power spectrum, Ce [4], or cquivalcntly 
the temperature angular correlation function, C(9). 

In Fig. 2 we also indicate specific curves which we se- 
lected based on the presence of particularly low variance 
annuli. This demonstrates that the simulated data con- 
tain similar low-variance circles as occur in the real data 
(e.g. the green curve could be claimed to have a 'family' 
of multiple low variance rings). The red curve was also 
chosen on the basis that it is systematically lower than 
average, similar to the GP10 I = 105°, b = 37° direction. 
This example was found by searching 500 random direc- 
tions, as opposed to 10885 in GP10. Hence it is clear 
that there is nothing unusual about this GP10 direction, 
and that it is simply an outlier with lower than average 
variance. 

Results for a simulated sky are also presented in GP10 
(the third panel of their figure 2), where the fluctua- 
tions in the variance appear very small compared with 
the real data. Details of this simulation are not provided 
in GP10. However, we found that we could reproduce 
a similar plot by using purely white noise simulations - 
i.e. Cg = constant, normalized to have approximately the 
correct total power. It is clear that the scatter around 
the average variance level for the GP10 simulation is not 
what is expected in a map containing the observed CMB 
anisotropics; instead, it agrees with that for a spectrum 
with entirely the wrong scale dependence. 

It is important to point out that even when a low- 
variance annulus is identified in real or simulated data, 
that does not imply the presence of a circular feature in 
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FIG. 2: Same as Fig. 1, except that this plot is for simulated 
W-band data. We have plotted the results for 100 random 
directions. The green and blue (heavy) curves were chosen 
from the full set of 500 random centres on the basis that they 
have annuli which show prominent dips in variance. The red 
(heavy) curve was also selected as it is systematically lower 
than average. It is clear that there are directions around 
which there are low variance circles, and even families of con- 
centric circles, just like in the real data. Again, high variance 
annuli could just as easily be highlighted. 



the map. The low variance could simply be due to par- 
ticularly low fluctuations in a segment of the annulus. 
Indeed, it is straightforward to repeat the above anal- 
ysis, but searching for low- variance features of shapes 
other than circular annuli. As an example, we performed 
the identical analysis steps, but looking instead for low- 
variance 'triangular annuli', i.e. the regions between con- 
centric equilateral triangles of different sizes [5]. The 
results arc presented in Fig. 3. It is clear from this figure 
that there are directions around which there are similarly 
low-variance triangles to the low-variance circles found 
above. Therefore there is nothing special about the pres- 
ence of low- variance circles on the sky. We also tested 
the same triangles in simulated data, finding results con- 
sistent with those in Fig. 3. 

We have refrained from any attempt to estimate the 
significance level for the specific patterns identified in 
GP10. Since the predictions for low variance annuli are 
quite vague, it would be impossible to avoid the issue of a 
posteriori statistics. But in any case, this is unnecessary, 
since it is clear that such low variance patterns occur by 
chance in realistically simulated skies. 

Recalling that the CMB 'signal' behaves statistically 
like a random 'noise', and that it is Gaussian in nature, 
it is possible to understand theoretically the behaviour 
of the variances in the annuli. All that is required is 
the Ce power spectrum, or its real-space equivalent, the 
2-point angular correlation function, which is essentially 
the variance on the sky as a function of angular separa- 
tion between positions. The variance of the temperature 
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FIG. 3: Same as Fig. 1, except that here we have calculated 
the variances in equilateral 'triangular annuli' inscribed into 
circles of radius R. We have plotted the results for 100 random 
directions and random orientations of the triangles. The green 
(heavy) curve was chosen from the random centres on the 
basis that it has several prominent low- variance triangles. 



Ti within an annulus, where i labels the pixel, is 

4 = ^E^ 2 -4E^ a) 



where N is the number of pixels. This is the variance for 
one particular realization of the Gaussian random field. 
Taking the expectation over all realizations in the ensem- 
ble (denoted (• • ■}) gives 
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In the limit of large R, we have (4) — > C(0), i.e. the 
expected variance is independent of the shape of the re- 
gion, and is given simply by the zero-lag angular correla- 
tion function C(0) = J2i( 2i + 1)Q/(4tt). This is around 
[110 A*K] 2 , which we estimate from the actual WMAP W- 
band map, which includes noise. This matches the levels 
seen in Figs. 1-3 for large R. 

Now we want the variance of o~\ over the ensemble. So 
we want 
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Straightforward computations give 
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FIG. 4: Theoretical values for the expectation (and ±ler lim- 
its) of the circular annulus variance (solid curves) compared 
to the mean and standard deviation from 1000 randomly cho- 
sen centres (dashed curves). To compute these we used a 
noiseless full-sky map simulated from the best-fit WMAP 7- 
year power spectrum [4]. Note that in the previous figures 
we show the standard deviation of the temperature in each 
annulus, as opposed to the variance. 



In Fig. 4 we show the theoretical values for the expecta- 
tion and variance of the annulus variance in Eqs. (2) and 
(5), compared to simulations based on a noiseless full-sky 
map from the best-fit WMAP 7-year power spectrum [4], 
convolved with a 13.2' Gaussian beam (which is close to 
the W-band beam size). We use only the CMB com- 
ponent for simplicity, since the ensemble average power 
spectrum (required to compute the theoretical values) is 
known in advance (whereas the actual W-band map also 
contains noise and a foreground component). This is why 
c 2 (R) in Fig. 4 is slightly below the levels in the previous 
figures. Note that Figs. 1-3 plot the standard deviation, 
a. to enable direct comparison with GP10, while Fig. 4 
plots the variance, a 2 . This is because the above deriva- 
tion does not carry over straightforwardly to the standard 
deviation. 

We see that the theoretical values agree very well with 
the results of actual calculations of the variance for many 
centres. The suppression in variance observed for small 
R is simply due to the presence of structures in the CMB 
on the prefered scale of the first acoustic peak at ~1°: 
the smallest annuli are roughly this size, and hence the 
temperature is expected to be more strongly correlated, 
i.e. the variances lower, than for annuli with larger R. 

In summary, the CMB probes the largest accessible 
scales, and so it is certainly worth looking for signatures 
which may be smoking guns for new physics (from the ex- 
tensive literature on this topic, we randomly suggest [6]). 
Such signatures certainly could reveal information about 
the earliest stages of the Universe. We have confirmed 
the presence of the low- variance circles in the CMB sky 
found by GP10. However, we have shown that such cir- 
cles are to be expected in typical realizations of the CMB 
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sky based on the usual assumption of Gaussian random 
primordial fluctuations with the scale dependence which 
is now so well measured. So if there are signals of ex- 
traordinarily early times buried in the CMB, they have 
not yet been found, and we will have to keep looking. 
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